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ABSTRACT 

A deep Chandra observation of the coohng core cluster Abeh 4059 (A4059) is presented. Previous 
studies have found two X-ray cavities in the central regions of A4059 together with a ridge of X- 
ray emission 20kpc south-west of the cluster center. These features are clearly related to the radio 
galaxy PKS 2354-35 which resides in the cD galaxy. Our new data confirm these previous findings 
and strengthen previous suggestions that the south-western ridge is colder and denser than, but 
in approximate pressure equilibrium with, the surrounding ICM atmosphere. In addition, we find 
evidence for a weak shock that wraps around the north and east sides of the cavity structure. Our 
data allow us to map the 2-dimensional distribution of metals in the ICM of A 4059 for the first 
time. We find that the SW ridge possesses an anomalously high (super-solar) metalicity. The unusual 
morphology, temperature structure and metal distribution all point to significant asymmetry in the 
ICM atmosphere prior to the onset of radio-galaxy activity. Motivated by the very high metalicity of 
the SW ridge, we hypothesize that the ICM asymmetry was caused by the extremely rapid stripping 
of metal enriched gas from a starburst galaxy that plunged through the core of A4059. Furthermore, 
we suggest that the onset of powerful radio-galaxy activity in the cD galaxy may have been initiated 
by this starburst /stripping event, either via the tidal-shocking of cold gas native to the cD galaxy, or 
the accretion of cold gas that had been stripped from the starburst galaxy. 

Subject headings: galaxies: abundances — galaxies: clusters: individual (A4059) — galaxies: individ- 
ual (PKS 2354-350) — galaxies: jets — X-rays: galaxies: clusters 



1. INTRODUCTION 

In the core regions of galaxy clusters. X-ray observa- 
tions directly reveal the radiative cooling of the intra- 
cluster medium (ICM). In most relaxed clusters, the ra- 
diative cooling time is measured t o be subst antially less 
than the Hubble time (see Peter son fc FabianI (|200 6)). 
In the absence of a balancing heat source, this cooling 
would result in a slow fiow of matter inwards towards 
the cente r of the cluster , a cooling fiow (|CQwie fc Binneyl 
[l977; Fab ian fc Nulsenlll977l : rFabianlll994f ). and a grow- 
ing central mass of cooled gas. The absence of signa- 
tures of this cold gas at non- X-ray wavelengths leads to 
the cooling fiow problem. More recently, high-resolution 
X-ray spectroscopy with XMM-Newton has shown an al- 
most complete lack of gas below a temperature of Tvir/3, 
where Tvir is the virial temperat ure of the c ore regions of 
the cluster ( Peterson et al. 200l: [Tamura et al. 2001). It 
is possible that gas below Tvir/3 ^1 — 2 keV is present 
but not detectable. For example, Fabian et al. (2001) 
discusses a model in which there are strong metalicity 
inhomogeneities and the high metalicity gas cools below 
1-2 keV too rapidly for detection. No compelling model 
for the creation of such inhomogeneities exists, however. 
Additionally, a genuine lack of cool gas would help ex- 
plain (and indeed may be required to explain) the defi- 
ciency of massive galaxies that is revealed while compar- 
ing the galaxy mass function and the dark matter halo 
mass function ([Benson et al.ll2003l ). 

The cooling fiow problem can clearly be solved if 
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some phenomenon heats the ICM core enough to, 
on average, balance radiative cooling. A leading 
possibility is that jets from active galactic nuclei 
(AGN) at the center of galaxy clusters somehow 
provide this heating. It has been well established 
that AGN undergo complex interactio ns with the 
surro u nding: I CM, resulting in c avities (iFabian et al.l 
20001: [Mc Namara et al. 2000; Blant on et all l200ll : 
Young et al. 2 002: Nulsen e t al. 2005b), ghost cavities 



(IFabian et alll2 000: McNa mara et all 120011: 1 Heinz et al 



'2002*: IChoi^eT al. 200 4), ripples (IFabian et all l200^ 
2005), shells ([F abian efap l2000t ). shocks (iJones et al l 
2004; Nulsen et al.ll2005al lbl). and filaments (Nulsen et all 
2005b). It is feasible that one or several of these inter- 
actions result in the required heating. 

In this paper, we present a new and moderately deep 
Chandra observation of the cooling core cluster Abell 
4059 (A4059; z = 0.049). Previous X-ray studies 
have shown that A4059 has a double peaked surface 
brightness distribution, with a bright ridge of X-ray 
emi ssion located to the southwest of the cluster cen- 
ter (iHuang fc Sarazinlll998l : iHeinz et al.ll2002l : IChoi et all 
12004). This ridge seems to be in pressure balance with 
the ambient material. At the western boundary of this 
ridge there is a sharp gradient in temperature and den- 
sity. It was suggested that this gradient could be ex- 
plained by radiative cooling resulting from an interac- 
tion between the ICM and the radio-galaxy induced dis- 
turbance. While the metalicity profile of A4059 was not 
clear from these previous studies, there was weak evi- 
dence that it peaks at roughly s olar values in the center 
regions and decreases outwards (|Choi et al.ll2004l ). 

Additionally, previous studies have shown two X-ray 
cavities in the center of A4059. The northwest cavity 
is clearly visible and distinct, while the southeast cav- 
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ity is more difEc ult to discern (|Huang fc Sarazi ^ [19981 : 
iHeinz et al.ll2QQ2f ). These previous studies found no evi- 
dence of shocks around the cavities. It was noted however 
that the axis connecting the centers of the two cavities 
does not pass through the position of the AGN itself, 
instead passing to the no rtheast of th e galac tic nucleus 
by - 10''. This prompted iHeinz et all (|2QQ2r ) to suggest 
that the radio-galaxy outburst erupted into an ICM that 
possessed a pre-existing bulk flow, and that the cavities 
have been pushed in the NE direction by that flow. 

High-frequency radio studies (1.4 — 8.5 GHz) of 
this cluster indicate the presence of two radio lobes 
(|Choi et al.ll2QQ4l : iTavlor et anil994D associated with the 
radio galaxy PKS2354-35. The northwest lobe only 
partially covers the northwest X-ray cavity, and the 
southeast lobe simply misses the southeast X-ray cavity. 
Both X-ray cavities are therefore ghost cavities, proba- 
bly formed from a previous powerful burst o f activity by 
the central FR I radio galaxy, PKS 2354-35 (|Heinz et al.l 
l2QQ2l : lChoi et al.ll2QQ4[ ). This has been recently confirmed 
by the detection of low frequency (330 MHz) radio emis- 
sion coincident with these cavities (T. Clarke, private 
communication) . 

As we discuss in this paper, our new Chandra data 
change our view of this cluster. Section [2] describes the 
Chandra observations and basic data reduction. Sec- 
tion [3] discusses additional data reduction and presents 
X-ray images, temperature and metalicity maps, and 
the results of the deprojection analysis. The results 
are discussed and interpreted in § [H In particular, 
we discuss how considerations of the metallicity asym- 
metry in A4059 leads us to a rather novel view of 
radio- mode feedback in this cluster. Finally, we present 
our conclusions in Section 5. Throughout this pa- 
per, we use J2000 coordinates and assume a standard 
WMAP cosmology (flat universe with Vt\ = 0.73, Hq = 
71kms-^Mpc~^; Spergel et al. 2003). Given the red- 
shift of A4059/PKS2354-35 of z = 0.04905, this gives a 
luminosity distance of 215.1 Mpc and a scale of 0.95 kpc 
per arcsecond (jWrightl 120061) ^. 

2. OBSERVATIONS AND DATA REDUCTION 

Chandra observed Abell 4059 on three separate oc- 
casions; 24-September-2000 (total exposure 24.6 ksec), 
4-January-2001 (total exposure 20.1 ksec), and 26/27- 
January-2005 (total exposure 93.6 ksec). Results from 
the first t wo, shorter, o bservations have alread y been pre- 
sented bv IHeinz et all (|2QQ2[ ) and lChoi etHI '(|2004f ). Ah 
observations were performed with the Advanced CCD 
Imaging Spectrometer (ACIS) placing the core of the 
cluster close to the telescope aim-point on the central 
back-illuminated chip (i.e., chip S3). All data were re- 
duced with CIAO version-3.4. The level-1 data were pro- 
cessed in order to apply the latest gain correction map 
as well as the latest version of the charge transfer ineffi- 
ciency (CTI) correction. The processed level-1 data were 
then filtered with a standard grade selection (including 
ASCA grades 0,2,3,4,6) resulting in new level-2 events 
lists. Periods of high background were searched for by 
examining the chip-Sl lightcurve. The September-2000 
data were found to be affected by a large background 
flare; elimination of the affected data reduced the good 

^ http:/ /www. astro. ucla.edu/~wright/CosmoCalc. html 



exposure time to 16 ksec. The 2001 and 2005 data were 
not affected by significant background flares, allowing the 
full 24.6 ksec and 93.6 ksec, respectively, to be utilized. 

The three datasets (totaling 134.2 ksec of usable data) 
were merged using the merge.all script in CIAO. A 4.2' 
X 4.2' (238 kpc x 238 kpc) region enclosing the clus- 
ter core was extracted to create an exposure-corrected 
image. It would, however, be inappropriate to use the 
merged dataset for spectral studies. The spectral re- 
sponse of the ACIS depends on the chip coordinates, 
which map to different physical coordinates for the three 
observations. Additionally, the spectral response has de- 
graded with time due to a slow increase in charge transfer 
inefficiency (CTI) and the build-up of a contaminant on 
the ACIS filter. For these reasons, spectral studies are 
only performed on the 2005 data. Given that the 2005 
observation corresponds to 70% of our total good data for 
this object, this restriction does not significantly impact 
the signal-to-noise of our spectral study. 

3. ANALYSIS AND RESULTS 

3.1. X-ray Images 

Figure [1] shows the exposure-corrected full-band (0.3- 
lOkeV) unsmoothed image of A4059. The double peak 
(corresponding to the cluster core and the SW ridge) is 
prominent, and the northwest cavity is clearly visible. A 
spur of bright emission is seen to the southwest of the 
core. The improved signal-to- noise allow previously un- 
known aspects of the morphology to be discerned. The 
core of the cluster and the SW ridge are clearly con- 
nected in the north by a high-contrast curved rim. This 
rim appears to continue through the brightest part of the 
SW ridge and define the southern parts of the ridge. Al- 
though much less distinct, there may be a similar struc- 
ture which is the approximate reflection of this rim in the 
line joining the cluster core and the brightest part of the 
SW ridge. At a much lower contrast, one can trace a thin 
plume extending from the core directly northwards for 
60 kpc (larcmin), and another plume extending north- 
west from the SW ridge. 

From the merged data, we produced adaptively 
smoothed, exposure-corrected, full band (0.3-10 keV) im- 
ages using the CIAO tool csmooth. Significance levels in 
the range 3 — 5<j were used to set the smoothing scales. 
The adaptively-smoothed image shown in Fig. [2] supports 
earlier morphological studies of Abell 4059. While the 
outer portion of the cluster is fairly relaxed, the inner 60 
kpc is rather disturbed showing the cavities and south- 
west ridge that have been well studied in previous works. 
Our deeper image reveals new structure, however. The 
adaptively smoothed image hints at part of a shell bound- 
ing the disturbed region of the cluster, especially to the 
southeast of the cluster. 

To investigate this apparent shell in more detail, we 
created an unsharp-mask image following the approach of 
Fabian et al. (2003). In detail, images with two different 
Gaussian smoothings were created from the merged data, 
one image with a Gaussian smoothing scale of 25 pixels 
(12. 5^') and the other with a Gaussian smoothing scale of 
3 pixels (1.5^'). The unsharp-mask image was created by 
subtracting the heavily smoothed image from the lightly 
smoothed image, a process which allows one to search 
for subtle discontinuities in either the surface brightness 
or the gradient of the surface brightness. The result is 
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Fig. 1. — Exposure-corrected full-band (0.3-10 keV) image of A4059. 

shown in Fig. [3l and reveals a subtle "edge" which has 
an approximate semi-circular form in the eastern and 
northern quadrants. The projected radius of this feature 
is 70-80 kpc. There are hints that, once this edge crosses 
the jet axis, this feature may sweep inwards and connect 
onto each end of the SW ridge. 

Another view of this structure is obtained by construct- 
ing radial surface brightness profiles. Figure [H shows the 
surface brightness profile for a sector of the cluster be- 
tween position angles 225° and 270°. This sector cuts 
the edge noted in the unsharp-mask image at a distance 
of 80 kpc from the cluster center and, indeed, a small 
discontinuity in gradient of the surface brightness can 
be seen at r = 80 kpc. Unfortunately, further studies of 
this feature (e.g., searching for temperature jumps asso- 
ciated with the surface brightness jump) are beyond what 
is possible with the current dataset. Rigorously charac- 
terizing the nature of this feature must await sigificantly 
deeper X-ray data. However, by analogy with findings 
in other systems (e.g., Virgo; Forman et al. 2007), it 
seems likely that this feature is a weak shock or com- 
pression wave driven into the ICM by the radio galaxy 
activity. In § 31 we discuss the implications of this hy- 
pothesis for the age of the AGN. As noted in Markevitch 



et al. (2000), a curved shock/compression front viewed 
in projection does not produce a true discontinuity in 
the surface brightness but, instead, is revealed through 
a discontinuity in the gradient of the surface brightness 
of the kind suggested by Fig. [H 

3.2. Imaging Spectroscopy Using Adaptive Binning 

By examining the spectrum of the X-ray emission 
across the image, the spatial distribution of temperature, 
metalicity, density, pressure, and entropy can be probed. 
This kind of physical information is crucial if we are to 
disentangle this complex source. In particular, under- 
standing the spatial distribution of metals in this cluster 
is a major motivation for these new data. 

Spectral analysis was performed on data that had 
been adaptively binned using the Weig hted Voronoi Tes- 
selation (WVT) binning algorithm of iDiehl fc Statler I 
(|2006l )^. This algorithm produces polygonal bins with 
a constant signal to noise ratio (SNR), where the num- 
ber of photons in each bin is the square of the prescribed 
SNR. Once the optimal polygon-tiling of the image plane 
has been determined, the CIAO tool specextract was 

^ littp://www.phy.oliiou.edu/~ diehl/WVT 
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Fig. 2. — Adaptively-smoothed, exposure corrected image of the core regions of A4059. Logarithmically-spaced contours of this same 
map are also over- laid. 



used to extract spectra, energy response matrices and 
effective area curves for each tile. A single background 
spectrum was obtained from the source free region at the 
north-east corner of the S3-chip. After experimentation, 
we determined that a SNR of 30 allows us to construct 
adequate maps of plasma temperature and emission mea- 
sure (and hence density, pressure and entropy). Metalic- 
ity maps requires higher quality spectra and could only 
be constructed once the SNR was increased to 40. 

The spectrum from each bin was modeled usin^ a single 
temperature thermal plasma model mekal (iMewe et al.l 
119851 . 119861 : iKaastral [19921 : [Liedahl et aH^lddb) modified 
by the effects of Galactic photoelectric absorption with 
column density A^h = 1-45 x lO^^cm"^. The abundance 
of metals within the plasma was a free parameter of the 
spectral models. However, the relative abundance of all 
elements heavier than helium was fixed at their cosmic 
ratios (Anders & Grevesse 1989); no statistical improve- 
ment in the fit is obtained by relaxing this assumption. 
Spectral fitting was performing using XSPEC version 
11.3.2. Although we formally fit the 0.5-10 keV data, 
the spectra of most tiles become background dominated 
about 7keV. The absorption column was held fixed for 
the spectral fitting presented here, although we have con- 
firmed that our results are robust to allowing Nn to be 



a free parameter. 

The temperature map (made from the SNR = 30 tes- 
selation) is shown in Fig. [5] (top panels). The most obvi- 
ous feature in this map is the region of cool (< 2.5 keV) 
gas that extends from the core through the SW ridge. 
Examination of this temperature map shows that cool 
gas has a close spatial correspondence with the highest 
surface brightness regions of the SW ridge. In particu- 
lar, the boundary of the region containing the cool gas 
closely follows the high-contrast curved rim of emission 
that connects the core with the SW ridge (see discussion 
in § 13. ip , and continues to track this curved rim through 
the SW ridge and into the southern extremity of the SW 
ridge. Interestingly, the spur that projects from the SW 
ridge in the north-westerly direction is appreciably hot- 
ter than the rest of the SW ridge (3.5-4.0 keV as opposed 
to 2.0-2.5 keV). 

For the first time, these data allow us to examine the 2- 
dimensional metalicity distribution for this cluster. The 
metalicity map derived from single temperature thermal 
plasma fits (using the SNR = 40 map) is shown in Fig. [5] 
(middle panels). The most significant feature in this 
map is the dramatic enhancement in metalicity (upto 
twice the average cosmic abundance) along the SW ridge. 
Metal enriched gas appears to be present along the en- 
tire length of the SW ridge, including the north-westerly 



5 




Fig. 3. — Unsharp mask image of A4059, created by subtracting a full band map smoothed with a cr = 12. 5^' Gaussian from a map 
smoothed with a cr = 1.5'' Gaussian. Overlaid are contours of the adaptively-smoothed map. 



spur which was noted above to be at a significantly higher 
temperature than the rest of the SW ridge. It is partic- 
ularly striking in this map that the metalicity of the SW 
ridge appears to be appreciably higher than that of the 
central core itself. 

The picture concerning the metalicity distribution in 
the centr almost regions changes somewhat, however, 
when we go beyond a single temperature thermal plasma 
model for the spectrum of each bin. In particular, we 
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Fig. 4. — Surface brightness profile for a sector of the cluster 
between position angles 225° and 270° (i.e., the ESE sector). 



have fitted spectra from the tiles of the SNR = 40 tes- 
selation with a two-component thermal plasma model; 
the temperature and emission measure of each plasma 
component is a free parameter of the fit, but the two 
components are assumed to share a common metalic- 
ity. Comparing the goodness of fit parameter for the 
single-component and two-component plasma models, we 
find that the two-component model is not a statistically 
better description of the data than the single-component 
model (employing the F-test for two additional model 
parameters) for the spectrum from most of the tiles, in- 
cluding those that make up the SW ridge. Furthermore, 
for this vast majority of tiles, the best-fitting value of 
the metalicity parameter is similar between the single- 
temperature and two- temperature fits. This last result 
is unsurprising since, in most cases, the normalization 
(emission measure) of one of the plasma components in 
these two-component fits vanishes, leaving essentially a 
single temperature fit. The spectra from the tiles in the 
centr almost bins of the cluster are, however, described 
significantly better (at the 90% or more level) by the 
two-component plasma model. The need for these multi- 
temperature models may reflect the morphological com- 
plexity of the ICM center or the combined contributions 
from both the core of the ICM and the ISM of the cD 
galaxy itself. 

It is important to assess the effect of the multi- 
component model on the inferred temperature and met- 
alicity distributions. We find that the one-component 
temperature maps shown in Fig. [5] (top panels) re- 
main an accurate approximation to the emission- 
measure weighted temperature map derived from the two- 




Fig. 5. — Top left panel : Temperature map of the core regions of A4059 based on the SNR = 30 Voronoi Tesselation. Color-bar 
units are keV. Top right panel : Smoothed temperature map, using a 20" box-car smoothing kernel. Middle left panel : Metalicity 
map (based on single temperature fits) of the core regions of A4059 based on the SNR = 40 Voronoi Tesselation. Color-bar units 
are cosmic abundances (as defined by Anders & Grevesse 1989). Middle right panel : Smoothed (single temperature) metalicity 
map, using a 40^^ box-car smoothing kernel. Bottom left panel : Metalicity map (based on two-temperature fits) of the core 
regions of A4059 based on the SNR = 40 Voronoi Tesselation. Color-bar units are cosmic abundances (as defined by Anders & 
Grevesse 1989). Bottom right panel : Smoothed (two-temperature) metalicity map, using a 40^^ box-car smoothing kernel. 
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Fig. 6. — The 0.5-10 keV spectra of an elliptical region en- 
compassing the SW ridge (upper spectrum; black points) and an 
off-ridge region of the same size located 12'' to the south-west of 
the SW ridge (lower spectrum; red points). Note that the iron-L 
complex at 1 keV and the iron-K line at ~ 6 keV are appreciable 
stronger in the ridge spectrum. 



component fit. On the other hand, the metahcity of the 
core inferred from the two-component fits is appreciably 
higher than that inferred from the single-temperature fits 
{Z = 1.2 - 2Zq compared with Z = 0.6 - I.OZq). The 
metalicity map derived from the two-component plasma 
fits Fig. [5] (bottom panels) still show a notable enhance- 
ment of metalicity along the SW ridge, but the metalicity 
enhancement now also extends into the regions immedi- 
ately surrounding the cD galaxy. 

To demonstrate the robustness of the conclusion re- 
garding the high abundance of the SW ridge, we have 
extracted the 0.5-10 keV spectrum from a single elliptical 
region (25'' x 10'') that encompasses the whole SW ridge 
and compared it with the spectrum from the same sized 
region located 12" to the south-west (i.e., beyond the 
SW ridge). One and two temperature thermal plasma 
fits to the SW-ridge spectrum implies a metallicity of 
Z = 1.70 ± O.I8Z0 and Z = IMtoir^e, respectively. 
On the other hand, one and two temperature fits to 
the non-ridge spectrum yields Z = 0.87 ± 0.24^© and 
Z = O.831q"26^0, respectively. Indeed, this metallicity 
difference is visually apparent in the spectra; examina- 
tion of the spectra (Fig. [6]) reveals an iron-L line complex 
at ~lkeV and iron-K line at ~6keV which are clearly 
stronger in the ridge spectrum than the off-ridge spec- 
trum. 

With some assumptions about the line-of-sight geom- 
etry, the emission measure can be used to compute the 
characteristic density of the plasma contributing to the 
observed radiation from a given spatial bin. Opera- 
tionally, the electron number density Ue is obtained from 
the normalization of the mekal component (N) in the 
one-temperature fits by 



Da{1 + z)\/4 X IOI^ttTVV-i, 



(1) 



where Da = 195 Mpc is the angular size distance to 
A4059, z = 0.049 is the redshift, and V is the effec- 
tive volume of the emitting plasma in that bin. The 
volume V was estimated by multiplying the area of the 
Voronoi tile by the distance to the cluster center, as de- 



fined by the location of peak brightness (which coincides 
with PKS 2354-35 itself). Given the density and spec- 
trally determined temperature, the pressure is obtained 
from the ideal gas equation: 

P = nkT. (2) 

where n is the total particle number density of the 
plasma, assumed to be given by n = 1.93ne (appropri- 
ate for 75% H and 25% He). The entropy index was 
calculated using 

5 = -pi- (3) 

where we assume a standard adiabatic index of 7 = 5/3. 
Finally, we compute the radiative cooling timescale 



'^cool 



3nfcT 



2nenjjA(T) 



(4) 



where uh is the hydrogen number density, and A(T) is 
the cooling function modeled using the fitting formula of 
Sutherland & Dopita (1993). 

These maps of density, pressure, entropy index and 
cooling time are shown in Fig. [3 They show the ex- 
pected trends; both density and pressure decrease with 
radius, while entropy and cooling time increases. It is 
particularly interesting to relate these physical quanti- 
ties to the SW ridge. Given the improve ment in the 
spatia l resolution of our maps over those of IChoi et al.l 
(|2004f ), we can spatially resolve variations of these phys- 
ical quantities ove r this ridge . Conf irming and strength- 
ening the result of Choi et al.l ((2004) , the pressure map is 
reasonably symmetric and does not show any discontinu- 
ity coincident with the SW ridge. Instead, the SW ridge 
corresponds to a region of denser, colder gas that main- 
tains a sharp interface (with a thickness of < 5 kpc) with 
the surrounding hotter gas. The gas within the ridge has 
a cooling time of 500-700 Myr, compared with cooling 
times in excess of 1.5 Gyr in the ICM it interfaces. 

3.3. Deprojection Analysis 

The 2-d spectral imaging analysis described above does 
not properly take into account the full 3-d nature of the 
cluster. Each bin will inevitably include contributions 
from gas that has a range of distances from the cluster 
center. Hence, any radial gradients of temperature, den- 
sity and metalicity will mean that each bin could there- 
fore include gases with different properties. 

To assess the effect of this, we performed a standard 
deprojection analysis using the pro j ct command within 
XSPEC. The observed cluster emission was divided into 
18 concentric, equal width annuli centered on the core 
of the cluster PKS2354-35. Assuming the cluster to be 
spherically symmetric (an assumption that clearly breaks 
down within the central region), we can "peel" spherical 
shells from the cluster; the spectrum from a particular 
shell has the contributions from more distant shells sub- 
tracted, thereby allowing one to model the spectrum of 
the gas only in that shell. We fitted each shell with a 
single temperature thermal plasma model (mekal) with 
absorption fixed at the Galactic value. The density of the 
gas in each shell was calculated from the emission mea- 
sure (which is directly proportional to the best-fitting 
normalization of the mekal model) by assuming that the 
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Fig. 7. — Maps of density (top-left), pressure (top-right), entropy index (bottom-left) and cooling time (bottom right). Color bar units 
are cm~^ (density), ergs"-*^ (pressure), Kcm~^ (entropy index), and years (cooling time). Contours of the adaptively smoothed full-band 
surface brightness are overlaid. 



plasma uniformly filled the shell. The pressure and en- 
tropy index follows straightforwardly once the tempera- 
ture and density are known. 

Th e results of this deprojection are shown in Fig. [HI 
As in 'Choi et al.l (2004), the temperature increases with 
radius, while the density and pressure decrease. The de- 
crease is is quite rapid along the inner radii. Between 10 
and 40 kpc, a slight dip in the pressure is visible, proba- 
bly due to the presence of the cavities. The deprojected 
metalicity results reflect the behavior seen in the spectral 
analysis of the 2-d (projected) data. The one-component 
plasma fit implies a metalicity that is low at the center, 
increases to a maximum at ~25kpc (corresponding to 
the metal rich SW ridge entering the annuli under con- 
sideration) and then decreases again at large radii. As 
in the case of the projected spectral analysis, this cen- 
tral "hole" in the metalicity is removed when a second 
plasma component (with a different temperature but in 
pressure equilibrium and with a common metalicity) is 
introduced (Fig. [9]). 

4. DISCUSSION AND CONCLUSIONS 

The jet axis of the central radio galaxy PKS 2354-35 is 
oriented on a SE to NW line, and this system has approx- 



imate reflection symmetry about a line passing through 
the AGN perpendicular to the jet axis; i.e., the jet and 
counter-jet sides look rather similar. However, this sys- 
tem shows striking asymmetries in the direction lateral to 
the jet axis. These new data uncover several new aspects 
of this asymmetry t hat were beyond t he r each of the 
previo us analyzes of iHeinz et"al] (|2002l ) and IChoi et all 

Hop). 

Firstly, both the unsharp mask image (Fig. [3|) and the 
radial surface brightness profile (Fig. [4]) reveals a sub- 
tle discontinuity in the gradient of the surface bright- 
ness that is approximately semi-circular in form (with 
60-80 kpc radius) and centered on the radio galaxy. The 
previous Chandra observations reported by (jHeinz et al.l 
2002; Choi et al. 2004) had insufficient signal to detect 
this feature. Although it cannot be rigorously proven 
from these data, it seems reasonable to interpret this 
as a weak shock driven into the previously undisturbed 
ICM by the onset of radio galaxy activity; such shocks are 
seen in all hydrodynamic simulations of jet/ICM interac- 
tions (e.g., Reynolds, Heinz & Begelman 2002; Vernaleo 
& Reynolds 2006, 2007 and references therein), and have 
been found by Chandra in Hydra A (Nulsen et al. 2005) 
and Virgo (Forman et al. 2007). Assuming that the ve- 
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radio-galaxy activity launched the shock, we have 
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Fig. 8. — Single temperature deprojection analysis for A4059. 

locity of this front through the ICM is greater than the 
local sound speed Cg and that the ICM is at rest with re- 
spect to the cD galaxy, we can estimate an upper limit on 
the time since the initial radio galaxy activity launched 
the shock. If L is the current distance between this shock 
front and the AGN, and t is the time since the onset of 
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Fig. 9. — Metalicity distribution obtained from the two- 
component thermal plasma model fitted to the deprojected data 
for A4059. 
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Using the deprojected temperature profile to determine 
Cs(r), we determine that t < 90Myr. This compares well 
with estimates for the source a ge based on the bu oyant 
rise time of the X-ray cavities (|Heinz et al.ll2QQ2l ). We 
note, however, that we only see this feature on the NE 
side of the jet line; there is no similar feature on the SW 
side. This may be due to the simple fact that the ICM 
emission on the SW side of the cluster is fainter than 
the NE side and, hence, the putative shock is below our 
detect ability limit. 

Secondly, the temperature map shows that the ICM 
temperature distribution is strongly asymmetric. A con- 
vex cold front (i.e., boundary across which there is a 
discontinuity in temperature but continuity in pressure) 
seems to extend from the central core through highest 
surface brightness portions of the SW ridge. This cooler 
and denser material has a comparatively short radiative 
cooling time, only 500-700 Myr. We note that this is 
rather lo nger than t he ce ntral/ridge cooling times in- 
ferred bv lChoi et all ([2OOI : we consider our cooling time 
map to be more robust given that we have the signal to 
noise to produce a map in which the g ross morphologica l 
structures are resolved [in contrast to IChoi et aH (|2004l ) 
who could only examine cooling times within bins that 
combined core and SW ridge emission]. One important 
qualitative difference is that our inferred cooling times 
are almost an order of magnitude longer than the plau- 
sible age of the radio source, eliminating the possibility 
discussed in lChoi et al.l (|2004l ) that the SW ridge is cool 
due to compressionally-enhanced radiative cooling. 

Thirdly, our new data allow us to map the 2- 
dimensional distribution of metals and reveal that it, too, 
is strongly asymmetric. High metalicity gas appears to 
extend along the length of the SW ridge, including into 
the NW spur of the SW ridge that lies beyond the cold 
front and may form the wall for the northern cavity. The 
result that the SW ridge possesses enhanced metalicity is 
robust in both one-component and two-component ther- 
mal plasma fits. We do find, however, that the apparent 
drop in metalicity within the central core that is inferred 
from one-component fits is removed when one considers 
two-component plasma fits. 

Assuming that the jet production process itself is ax- 
ially symmetric, all of the evidence points to an ini- 
tial (pre-interaction) ICM atmosphere that was strongly 
asymmetric. The asymmetries in this system have pre- 
viously been interpreted in terms of a bulk flow within 
the ICM (flowing in the plane of the sky from the SW 
to the NE) prior to the onset of powerful radio galaxy 
activity. However, it is very hard to understand the met- 
alicity asymmetry in terms of such a bulk flow. Instead, 
it seems likely that an anomalously dense, cool, metal- 
rich region of gas was situated to the SW of the center of 
the cD galaxy at the time that the AGN initiated its jet 
activity. This anomalous region would have been weakly 
shocked and possibly pushed higher into the ICM atmo- 
sphere by the laterally expanding radio galaxy cocoon, 
thereby forming the SW ridge. If this scenario is cor- 
rect, it seems that some fraction of this metal rich gas 
has been entrained into the walls of the X-ray cavities. 
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especially in the NW direction. 

At this point, it is useful to estimate the mass of gas 
within the SW ridge and the gravitational potential en- 
ergy of this structure. The brightest region of the (pro- 
jected) SW ridge emission can be encompassed by an el- 
lipse with major and minor axes of 15 arcsec and 8 arcsec 
respectively. Assuming that this is a oblate spheroid 
in projection, we estimate the emitting volume to be 
V = 2 X 10^^ cm^. A single-component thermal plasma fit 
provides an adequate description of the data and implies 
a total emission measure of 

EM = [ nldV ^Ix 10^^ cm-^ (6) 
Jv 

Assuming that twice cosmic abundance plasma uniformly 
fills this region, we conclude that the total mass of ICM 
within this ridge is M ~ 5 x 10^ Mq. We can also esti- 
mate the energy required to lift this matter from a loca- 
tion close to the center of the cD galaxy to its current 
position in the gravitational potential. Assuming that 
the undisturbed ICM is approximately isothermal with 
sound speed Cg and is in a hydrostatic configuration with 
density profile p(r), integration of the equation of hydro- 
static equilibrium gives the gravitational potential. 



where 7 = 5/3 is the ratio of specific heat capacities 
and po is some fiducial density that plays the role of 
the (uninteresting) constant of integration. The change 
in gravitational potential energy when mass M is lifted 
within the atmosphere is given by 

= MA$ = ^ In f (8) 

7 \P2j 

where pi and p2 are the densities of the undisturbed 
ICM atmosphere at the start point and end point, re- 
spectively, of the lifting process. Using values of Cs ~ 
700kms~^ (corresponding to to the kT ^ 2keV plasma 
that characterizes the central regions of the ICM), and 
Pi/p2 ^ 5 (from the deprojection results), we obtain 
AE ^ 5 X 1 0^^ er g. To place this energy into context, 
iHeinz et"aD (|2QQ2f ) estimated that the energy required to 
inflate the observed X-ray cavities exceeded 8 x 10^^ erg, 
almost 20 times greater than the energy required to lift 
the cool gas in the SW ridge. 

While it seems clear that the ICM atmosphere was 
asymmetric prior to the (most recent) onset of radio 
galaxy activity, the cause of that asymmetry is unclear. 
One possibility is that the asymmetry was caused by the 
merger /assimilation of A4059 with a small galaxy cluster 
or group. The SW ridge could then be identified with the 
remains of the ICM core of the smaller subcluster (albeit 
with a morphology that has been subsequently shaped by 
interaction with the jet-blown cocoon of PKS 2354-35). 
Indeed, Heinz et al. (2003) have shown that some frac- 
tion of ICM cold fronts are likely to be associated with 
subcluster merger events, and that the hydrodynamics 
of such an interaction can cause high metalicity gas to 
be dredged up from the center of the merging core and 
carried up to the cold front. This would qualitatively 
provide a natural explanation for the high metalicity of 
the SW ridge. The gross X-ray morphology of A4059 is 



also reminiscent of the "X-ray comet tail" found in the 
Abell 2670 cluster (Fujita, Sarazin & Sivakoff 2006). The 
mass of the gas in the "comet tail" is similar to the SW 
ridge in A4059, and it also appears cooler than the am- 
bient ICM. In the case of A2670 the association of the 
"comet tail" with a large elliptical galaxy, and the high 
peculiar motion of the cD galaxy, makes it very plausi- 
ble that we are witnessing the merger of a small cluster 
or group with A2670 and that the comet-tail is the dis- 
torted/stripped remnants of the subcluster 's ICM core. 

The subcluster merger scenario faces two potential 
problems. Firstly, unlike the case of A2670, there is 
no obvious association of the SW ridge with any galaxy. 
Since the galaxies associated with the merging subcluster 
may have left behind their ICM core (such as is seen in 
the Bullet Cluster IE 0657-56; Markevitch et al. 2002), 
this fact by itself does not rule out the merging subclus- 
ter hypothesis; it does predict, however, that a detailed 
mapping of the galaxy dynamics in this cluster should 
reveal a kinematically distinct subcluster. Secondly, the 
metalicity of the SW ridge is 1.5 — 2 Z©, an extreme value 
if this is to be identified with the ICM core of a merg- 
ing subcluster (Baumgartner et al. 2005; Snowden et 
al. 2007). Hence, while the merging subcluster hypoth- 
esis cannot yet be ruled out, we are led to explore other 
scenarios. 

The very high metalicity of the SW ridge, in particular, 
pushes us to explore scenarios in which the gas originates 
directly from the metal enriched ISM of an actively star- 
forming galaxy. In order to deposit 5x10^ M© of metal 
rich gas into a localized region of the cluster, the donor 
galaxy must have been massive (comparable to or greater 
than an L* galaxy) and the stripping mechanism must 
have been efficient and rapid. One possibility is that a 
large gas-rich late-type galaxy plunged through the core 
of A4059 for the first time and underwent an intense 
compressionally-driven starburst (Reverte et al. 2007). 
The combination of the starburst-driven superwind and 
ram-pressure stripping by the ICM could have allowed 
the galaxy to dump the required amount of metal-rich 
gas on a timescale short compared with the time taken 
to traverse the cluster core. Essentially, the star-bursting 
galaxy would leave a "debris" trail behind it compris- 
ing of metal rich, relatively cool, possibly multiphase 
gas. We hypothesize that this debris trail is the high- 
metalicity inhomogeneity which, once impacted by the 
radio-galaxy activity, produces the SW ridge feature in 
A4059. 

If the projected (i.e., plane of the sky) velocity of 
the galaxy during its passage through the cluster core 
is V = 1000^3 km s~^, the metal rich ISM would have 
to be stripped from the galaxy at a rate exceeding 
200i'3 M0yr~^ in order produce a distinct region consist- 
ing of 5 X 10^ M0 of metal-rich gas that is only ~ 25 kpc in 
projected extent. Once stripped away from the galaxy, 
this denser gas will follow a decaying orbit within the 
(dark matter dominated) gravitational potential due to 
the action of gas drag. Concurrently, hydrodynamic in- 
stabilities will tend to mix this stripped gas with the 
ambient ICM. Detailed theoretical modeling beyond the 
scope of this paper is required in order to estimate the 
time over which the starburst-delivered gas will remain 
a distinct inhomogeneity within the ICM core. For now, 
we estimate an upper limit on the time since this star- 
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Fig. 10. — Digitized Sky Survey image of a 10 x 10 arcmin^ field in the central regions of the A 4059 cluster. The cD galaxy ESO 349-GOlO 
hosts the radio-galaxy PKS 2354-35 and resides at the center of the ICM atmosphere. The bright spiral galaxy ESO 349-G009 can be 
clearly seen to the north of the cD galaxy. 



burst /stripping event occurred by equating it with the 
orbital time of the gas within the potential, r ~ ^-K^JVJg^ 
where g = — V<l> is the acceleration due to gravity in the 
cluster potential. This timescale is within a factor of or- 
der unity the same as the time taken for the cooler gas 
to fall radially into the center of the potential. Assuming 
that the underlying ambient ICM density field is approx- 
imately spherically symmetric p = p(r)^ we have 



and we see that 



r ~ 27r\ 



-rdp/pdr) ' 



(10) 



Putting r = 25kpc, and noting that —rdp/pdr ~ 1 — 2, 
we estimate that the starburst/stripping event occurred 
at most r - 200 - 300 Myr ago. 

An obvious and important prediction of this starburst- 
stripping hypothesis is that the responsible galaxy should 
still be in the vicinity and should still display post- 
starburst signatures. Using the above estimate of the 
upper limit on the time since the starburst/stripping 
event, the post-starburst galaxy should still be within a 
projected distance of SOOT'S kpc of the cluster core. This 
corresponds to an angular distance of 285'L'3 arcsec. Of 
course, deceleration of the galaxy as it climbs out of the 
dark matter potential would reduce this distance. 



A systematic study of the optical spectrum of the mem- 
ber galaxies of A4059 is required in order to find the 
post-starburst galaxy predicted by this hypothesis. One 
interesting candidate, however, is the bright spiral galaxy 
ESO 349-G009 situated 280'' directly to the north of the 
cD galaxy ESO 349-GOlO (see Fig. [TO]). The radial ve- 
locity difference between ESO 349-G009 and ESO 349- 
G010/PKS2354-35 is almost = 2100 km s"^ {v = 
12612 kms-^ for ESO 349-G009 compared with v = 
14705 kms-^ for ESO 349-G010/PKS2354-35; de Vau- 
couleurs et al., 1991) implying that this late- type galaxy 
is either simply a foreground object or a high-velocity 
galaxy close to the A4059 cluster. From the absolute 
2MASS K-band magnitude of Mk = —25.82, we can as- 
sume a K-band stellar- mass to light ratio of 0.7 (Bell et 
al. 2003) and a solar K-band absolute magnitude of 3.27 
(Holmberg et al. 2006) in order to derive an estimate 
for the stellar mass of this galaxy, ^ 3 x 10^^ M©. 
Thus, we see that the SW ridge in the ICM core could 
have been formed via the ejection of gas totaling a few 
percent of the stellar mass. 

What makes ESO 349-G009 a particular interesting 
candidate is the presence of X-ray and near ultravio- 
let (NUV) evidence for on-going vigorous star forma- 
tion. Sun et al. (2007) present an analysis of the X- 
ray emission from ESO 349-G009 using the same Chan- 
dra data as reported upon here. They find a lumi- 
nous (spatially resolved) X-ray corona associated with 
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the galactic bulge in addition to three point sources as- 
sociated with the spiral arms (Fig. [TT]). They also note 
that, despite the absence of an AGN (i.e., an unresolved 
source at the location of the galactic nucleus) in this 
galaxy, the X-ray emission has a hard component with 
^2-10 /ceV = 0.8 — 2.7 X 10^-^ ergs"-^. Associating the 
hard X-rays with High Mass X-ray Binaries, they use 
the correlation of Grimm, Gilfanov & Sunyaev (2003) to 
estimate a star formation rate of > 12 MQyr"^. This con- 
clusion is supported by the fact that ESO 349-G009 is lu- 
minous in the NUV band (see Fig. 26 of Sun et al. 2007). 
Within the context of our starburst-stripping hypothesis, 
one can postulate a scenario in which ESO 349-G009 ex- 
ecuted a high velocity plunge through the core of A4059 
and suffered an intense ICM-pressure induced starburst 
which is still in the process of dying away. High-quality 
optical spectroscopy of ESO 349-G009 combined with 
stellar-synthesis models will allow its recent star forma- 
tion history to be constrained and will provide a powerful 
check of this scenario. In particular, an intense starburst 
100-300 Myr ago will be revealed through strong Balmer 
line absorption features resulting from the large numbers 
of A stars. 

It must be noted, however, that even just the cur- 
rent line-of-sight velocity difference between ESO 349- 
G009 and ESO 349-G010/PKS2354-35 likely requires 
that ESO 349-G009 is marginally unbound to A 4059. 
Thus, if this galaxy is indeed the origin of the high metal- 
licity plume, this system is a rare example of an unbound 
massive galaxy falling through the core of a rich cluster. 

The temporal coincidence between the onset of 
powerful radio-galaxy activity and the hypothesized 
plunge/stripping of a large starburst galaxy suggests a 
connection. There are two principal possibilities. The 
close passage of a high velocity massive galaxy would 




tidally-shock any cold gas within the cD galaxy, driving 
an inflow and possibly triggering AGN activity. Alterna- 
tively, the gas stripped from the starburst galaxy is very 
likely to be multiphase (with regions of cold, atomic gas 
entrained by the hot X-ray emitting gas). Once stripped 
from the passing galaxy, the cold and dense atomic gas 
would preferentially fall (or, more precisely, sediment) 
into the cD galaxy and, again, possibly fuel the AGN ac- 
tivity. Either way, this presents a very different model of 
"radio-mode" AGN feedback to that normally envisaged; 
in this case, we are witnessing AGN feedback triggered 
directly by the close passage of a gas rich galaxy rather 
than radiative cooling of the host ICM. 



Fig. 11. — Digitized Sky Survey image of a 2 x 2arcmin^ field 
centered on ESO 349-G009. Overlaid are contours of the full-band 
(0.5-10 keV) X-ray intensity. See text and Sun et al. (2007) for a 
discussion of this emission. 
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